Mapping a depinning transition to polynuclear growth 
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We develop a phenomenological mapping between submonolayer polynuclear growth (PNG) and 
the interface dynamics at and below the depinning transition in the Kardar-Parisi-Zhang equation 
for a negative non-linearity A. This is possible since the phase transition is of first-order, with no 
diverging correlation length as the transition is approached from below. The morphology of the 
still-active and pinned configurations and the interface velocity are compared to the PNG picture. 
The interface mean height scales as erf(t). 
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I. INTRODUCTION 

The dynamics of driven manifolds in random media 
presents many examples of non-equilibrium phase tran- 
sitions. The interest hes often in how an object inter- 
acts with a quenched disorder environment. If the driv- 
ing force is small enough and the temperature zero, the 
manifold, e.g. a domain wall in a magnet, or an interface 
between two phases, can get pinned. This means, quite 
simply, that its velocity becomes zero. Therefore one can 
discuss the physics in terms of an order parameter (ve- 
locity) and a control parameter (external force). At and 
close to the critical force value Fc the driven interface 
can develop critical correlations due to the frustration 
that arises from the competition between, typically, elas- 
tic forces and the randomness of the medium [Q-^ . 

An example is provided by the quenched Kardar- 
Parisi-Zhang (QKPZ) equation, 

= '^y'h + ^{vhr +F+ r?(x, h). (1) 



where h is the interface height, a surface tension, A 
the strength of the nonlinearity characterizing the KPZ 
class 0], _F a driving force, and ij is the quenched noise. 
If A is positive and stays finite in the limit F F^ , 
the depinning phase transition in 1+1 dimensions can be 
mapped to directed percolation (DP) via a geometrical 
construction. In this way one can compute the critical 
exponents through the knowledge concerning the DP ex- 
ponents, and the particular model where this is explicitly 
possible has been coined 'directed percolation depinning' 
one for obvious reasons ||^~p^. 

If A is negative, however, the character of the transi- 
tion changes completely. Fig. 1 shows an example of an 
interface at F ^ Fc. The depinning in this case was in- 
vestigated by Jeong, Kahng, and Kim Q, who showed 
that the transition becomes of a first-order kind and that 
this persists even in the presence of finite tilts below a 
critical rric. The aim of this article is to use a new ge- 
ometrical construction to discuss the morphology in the 
case depicted in Fig. 1, similarly in spirit to the DPD- 
mapping of the A > 0-case. 

The essential physics in the A < QKPZ class, at and 
below Fr arises from "facet formation" , the fact that the 



slopes of the triangular mounds are just dependent on 
the parameters of the equation, and from "valley forma- 
tion". A look at Fig. 1 where the interface is shown at 
various times ti implies that it is these two mechanisms 
that are responsible for the final shape of an interface that 
gets pinned. In this article we map the two-dimensional 
growth of the QPKZ interfaces to a one-dimensional sub- 
monolayer polynuclear growth model (PNG) ||l^ based 
on this picture of a statistically constant growth velocity 
of the facets, and random pinning events on the roughly 
speaking planar areas of growth. This implies that one 
can understand the pinning of an interface in a quenched 
environment by considering a model with 'thermal noise', 
only. The reason why this is possible is the first-order 
character of the transition: there is no diverging corre- 
lation length, and therefore the physics of the depinning 
can be understood based on a local picture of events. No- 
tice that the PNG in the multilayer version is actually in 
the thermal KPZ-class, and that the relation in our case 
between a growth model and the PNG is indirect. 
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FIG. 1. Temporal evolution of a QKPZ interface at an ex- 
ternal driving of _F = 14.5. The snapshots have been taken 
every 1000th time step. The vertical slicing lines designate the 
positions of "valleys" and "peaks" located by the algorithm 
described in the text. 
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Section two of this paper introduces the actual numer- 
ical QKPZ model and the relevant parameters. In Sec- 
tion three, we discuss in detail the mapping between the 
QKPZ and the PNG models including a measurement of 
the model parameters for use in the PNG model. We 
also discuss the relevant parts of the PNG literature. In 
Section four the free parameters (growth velocity, pin- 
ning probability per unit time) are discussed and com- 
puted, and a basic comparison of the QKPZ interfaces 
with those that come from the PNG representation is 
presented. Section five finishes the paper with a discus- 
sion. 



II. NUMERICS 

The model we focus on is a 1 -I- 1 dimensional QKPZ 
system with the spatial coordinate x and time t dis- 
cretized, while the interface height /i is a continuous num- 
ber. We use periodic boundary conditions with flat ini- 
tial configuration. The forthcoming results have been 
obtained from systems in which the equation coefficients 
are invariably v = 10, \ — —4, and g ~ 20. 

The system is updated simultaneously for all lattice 
sites in every time step according to the finite differences 
version of the governing equation (|l|) with = 0.01: 

f,(t) = v[h,+i{t) + h,^i{t) - 2h,{t)] + 
+ - h,^i{t)f + F + gT^,{[h,{t)]) 

M^ + At) = |J1;! + ^'^^' ^fi'W.>0 (2) 
^ ' \hi{t) otherwise ^ ' 

The value of the noise term rji at integer positions of 
the height field denoted by [•] is chosen randomly from 
the interval [—1; 1] with a uniform distribution and is un- 
changed in time. The choice of the noise term is meant 
to mimic the simulation carried out in pT[ . Neverthe- 
less, the particular form of rji has a mere quantitative 
influence on some measurables such as the critical driv- 
ing force Fc] the qualitative behavior of the systems have 
been checked to yield consistent results independently of 
the various options for the representation of rji. 

We study the morphology of slightly subcritical inter- 
faces with the mapping to the PNG model in mind. To 
this end we identify the macroscopic facets along the in- 
terface locating the "valleys" and "peaks" in the entirely 
evolved and stopped systems. 

The algorithm to find the extrema is composed as fol- 
lows. As a first step the pinned interface is scanned for all 
of the local minima and maxima positions. Subsequently 
after that, the pair of adjacent minimum/maximum with 
the smallest distance to each other is removed and this 
step repeated until the distance between the two closest 
extremal positions is no more than a specific parameter. 
Applying this procedure we eliminate the noise- like small 
height fiuctuations of the interface. The resolution of the 



facets depends on the actual choice of the smallest dis- 
tance parameter and we use a few lattice units for that. 



III. MAPPING TO THE POLYNUCLEAR 
GROWTH MODEL 

A. The interface growth in one dimension 

Following the time evolution of some interface config- 
urations like the one in Fig. |l|, it suggests itself that the 
pinning events resulting in the valleys of the final mor- 
phology take place in a seemingly random way on the 
growing terraces. The interface is observed to propagate 
with a uniform average vertical velocity at the plateaus 
of the facets, in other words, the height of the separated 
terraces of the interface is on the average on the same 
level at all times. The straight sides of the facets mean 
that the horizontal distance between the end points of 
neighboring facets at the level of propagation grows in a 
linear fashion with time, following to first order from a 
balance between the nonlinearity and the driving force in 
Eq. (I). 

If we consider the h axis of Fig. |l| as time and the 
sides of the final interface facets as phase boundaries in 
this model, we can outline a picture of a model mapping 
as follows. We define a one dimensional lattice of cells 
on the grounds of the PNG model with two-state filling 
possibilities: a site is either occupied by a particle or 
empty, with an initially empty configuration. During the 
simultaneous updating steps of the lattice, a seed par- 
ticle may be deposited to any of the empty cells with a 
certain probability in every time step. The actual deposi- 
tion probability and place may depend on various factors 
which will be discussed later. 
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FIG. 2. Time evolution of PNG islands with the seeds filled 
black. The mapping is performed for the interface configura- 
tion depicted in Fig. ^ with a system size rescaling factor of 
A = 0.05. 
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Deposition of seeds creates growing islands with one 
particle added per time step at both ends according to 
the formal rule of 



1 if Si_i(i) = 1 or St+i{t) = 1 
Si (t) otherwise 



(3) 



Si{t) denotes the filling state of the cell i at time t. 
An occupied site is marked by 1 while an empty one is 
marked by 0. In other words (||) means that a cell gets 
filled if either of its neighbors is already filled. The sys- 
tem has periodic boundary conditions so Si-i — sl for 
i = 1 and s^+i — si for i = L. L is the lattice size and 
the indices run from 1 to L. 

It is obvious that if we associate the empty cells in 
this model with the still-growing parts of the interface 
and the filled cells with the uncovered regions, the time 
evolution of the system would give a very similar picture 
to the final configuration of Fig. |l|, if the way the seeds 
are deposited is chosen the right way. An example of the 
mapping is demonstrated by Fig. |[ 

Fig. H suggests that the probability the first seed is 
introduced at a given time should decay exponentially 
in time, which is the earmark of a geometric process. 
This means that the probability of a seed deposition 
(the creation of the first pinned site) has a constant 
probability per each time step. If this is true, then 
the growth closely resembles the submonolayer polynu- 
clear growth model (PNG), or the Kolmogorov-Avrami- 
Johnson-Mehl (KAJM) model (The PNG model 

is originally defined so that the growth itself takes place 
in many layers atop each other.) 
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FIG. 3. The probability distribution P of the first pinning 
height h in the KPZ system, rescaled by the system size of 
500, 1000, and 2500 for O, A, and □, respectively. The driv- 
ing force is F = 14.8. The averaging is done for 10000 inter- 
face configurations in each case. The rescaled distribution is 
found to decay exponentially with an associated exponent of 



Another crucial question is: how is the pinning prob- 
ability or the nucleation probability in the PNG picture 
dependent on F7 Fig. ^ illustrates the progress of two 
quantities close to the estimated Fc- Both the average 
number of pinning events and the average height (time) 
of the first pinning event scale roughly as power laws with 
\F — Fc\. Since the depinning transition is of first order, 
it is not a priori obvious how these quantities should be- 
have as F^ . The most natural guess is that the 
pinning probability per site per time-step goes to zero 
continuously, which would indicate that h{F) and N{F), 
as depicted in Fig. |[ should approach infinity and zero, 
respectively with cut-offs imposed by the maximum time 
the interfaces are followed, and the system size, respec- 
tively. 
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FIG. 4. The first pinning height h as function of the driv- 
ing force \F — Fc\/Fc for a system of 500 lattice units. The 
saturation is caused by limitations imposed on the maximum 
time configurations are allowed to evolve because of practi- 
cal reasons. The inset shows the number of pinnings taking 
place. 



B. Interface morphology 

The valleys and peaks are identified like described in 
Section II. The side gradients of the facets are measured 
on both sides and are displayed for various system con- 
figurations in Table | The height and position of the 
pinning events are also determined and recorded. F is 
chosen so that the system develops only a limited num- 
ber of facets; this way, the uncertainty of facet identi- 
fication and measurement stays under control. On the 
other hand, we are mainly interested in the behavior 
when F ^ Fc- We do not determine Fc with very great 
accuracy; good estimates can be found in Q . 

According to Table |, the facet side gradients are inde- 
pendent of system size for a given F . Additionally, the 
left and right hand side slopes are the same as well, as 
expected in the absence of a global tilt. In contrast to 
this, the number of pinning events taking place in the 
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systems is subject to finite size scaling. The error values 
associated with the averages are apparently large (30- 
40%). 

C. PNG theory 

For the next section where the growth statistics are 
compared with the one-dimensional picture some results 
are needed about the interesting quantities. In the ordi- 
nary KAJM theory islands are nucleated on a line ran- 
domly, with a rate 7 and grow at a constant velocity. 
To map to the KPZ case one needs a function that re- 
lates that velocity to the KPZ mound slopes as discussed 
below. If one has N{t) nucleated islands (here pinned 
parts of the interface) then the fraction of uncovered 
space (here growing fraction of interface) S{t) satisfies 
the rate equation dS/dt = —N with a proportionality 
factor chosen as unity. It is easy to understand that in 
a large system the temporal fluctuations of the slopes 
of the growing KPZ interface average out, to an average 
growth velocity in the KAJM picture. This holds if the 
slopes do not depend on e.g. the size of the terrace they 
are adjacent to which seems trivially true. Then the only 
factor that needs to be checked is whether the nucleation 
events happen randomly, with a constant density as in 
the simple KAJM as discussed above. 

In the one-dimensional picture, one can proceed by 
defining size distributions for the islands, for the gaps 
between the island, and finally the density of islands that 
contain a fixed number of seeds, n. The last one implies 
that such islands result from the coalescence of n single- 
seed islands. Such distributions describe exactly the KPZ 
interface in the KAJM language (assuming that the map- 
ping to PNG/KAJM works out). The last of these dis- 
tributions is hard to compute since it involves a solution 
of a Smoluchowski-like equation for the densities for dif- 
ferent n. We quote here some of the appropriate results 
from |l^Jl^. First, the island number density, 

N{t)r^te-'\ (4) 

corresponds to the number of valleys in the KPZ picture 
(as, also, a function of average interface height, N{h)). 
Second, the uncovered fraction gives directly the inter- 
face velocity after a proper normalization, and reads 

S{t) = (5) 
implying a fast decay of the growth. 

IV. COMPARISON OF THE PNG PICTURE 
WITH THE KPZ MODEL 

Following the considerations above, let us denote the 
probability of a seed deposition per lattice unit in the 
PNG model 11 for every time step. 11 should be tuned 



appropriately so that the average number of seed depo- 
sitions equals the average number of pinning events or 
valleys in the KPZ system. This is the only parameter 
that affects the "morphology" in this system, while an- 
other one, an artificial "side gradient" ^ has also to be 
introduced to transform the time scale of the PNG model 
to length scale in the KPZ counterpart. 

Notice that the virtual valleys and summits from 
the PNG model after the scale transformation are not 
rounded. The presence of the surface tension term in the 
KPZ equation counteracts such sharp corners, though. 
Still, if the strength of the surface tension coefficient ly 
is sufficiently small, the resulting interface is very similar 
to the one obtained from the PNG model. 

To achieve proper quantitative predictions and match 
between the KPZ morphology statistics and the PNG 
picture, we need to tune 11 and /j, so that measuring the 
first pinning height (seed deposition time) distribution, 
we end up with the same decay exponents. We relate the 
general height in the KPZ model h^pz to time in the 
PNG system tpj^G as Hkpz = ^-tpNG- The probability 
of the first seed deposition follows a geometric process: 

P{tpNG =t) = P{hKPz/l^ = r) = (1 - py-^p = 

= (1 -p)X'^-''^/''-V 

log P[tpNG =t)= log \ XKPZ , (6) 

l-p ^ 

if p denotes the seed depositing probability for the entire 
length of the PNG system, that is, p — HLpng- Xkpz is 
the first pinning height (the altitude of the lowest valley) 
in the corresponding KPZ system. Comparing Fig. || to 
Eq. (|^) we find that the assumption of p — const, is valid. 

The experimental value of the probability decay expo- 
nent /3 in the KPZ system should thus equal log(l ~p)/ jj. 
expressed by the PNG parameters. The Taylor series ex- 
pansion of the later results in p = — /3/ilnlO up to first 
order. 

If the horizontal scale transformation between the KPZ 
and PNG systems is given by LpjvG — ^Lkpz, then fi 
can simply be calculated from the measured facet side 
gradients in the KPZ system, summarized in Table |. 
Then if s denotes the empirical value of the side slope, 
fi = s/A. According to these, the normalized seed depo- 
sition probability 11 is obtained as 

/3s In 10 

ALpNG 

The discussion above assumed that the parameter 
matching is performed between two particular lattices 
of specific features, most importantly the number of cells 
in the systems. Considering the data collapse of Fig. ||, 
though, an additional relation for /3 emerges instantly: 
f3 ~ Lkpz, provided that the equation parameters (v, 
A, F, g) are otherwise fixed. Since s does not explicitly 
depend on Lkpz, and \i (3 — bLxpz, 



4 



n 



bLKPzslnlO 6s In 10 



A{ALkpz) 



A2 



(8) 



A further consequence of the fact that P ~ ^(hLxpz) 
for the first pinning probability with a suitable scaling 
function is that the general pinning probability per 
time unit in propagating fronts of a QKPZ system is lin- 
early proportional to the integrated local interface veloc- 
ity, v{x,t), 



p 

J n'. 



v{x, t) dx. 



(9) 



Using the parameters above it is possible to completely 
recover the exact scaling law of the first pinning event in 
a system conforming to the KPZ equation. It may still 
be uncertain, though, whether the subsequent pinning 
events follow the same governing rules as the first one, in 
particular the assumption that the seed creation proba- 
bility is linearly proportional to the island length. 
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FIG. 5. The "island number" distribution N{h) calculated 
for a KPZ system according to the text. Since the inter- 
face propagation velocity is constant, h can be appropriately 
rescaled so as to yield time in a PNG system. The system 
size is 2500, while the driving force is F = 14.5. The num- 
ber of averaging iterations is 10000. The theoretical func- 
tion of N{h) = Ahexp{—Bh^) is fitted to the empirical data 
points minimizing the quadratic error and giving the results 
oi A = 5.84 • 10"^ and B = 2.33 ■ 10"^ The slight shift of 
the distribution maximum may come from the fact that an 
initially flat interface cannot get pinned before the the nec- 
essary fluctuations develop. On the other hand, the omission 
of very small facets can have an unwanted effect on the data 
as well. 



the interface, reducing the actual number of islands by 
one. An illustration of the island construction is shown 
in Fig. ^. This thought experiment leaves us with a way 
to construct the associated island number distribution 
in function of the time (height) in the KPZ model. The 
N(t) function for the individual samples is naturally com- 
posed of step functions taking on the appropriate integer 
constants between the respective time (height) intervals. 
The PNG and KPZ behaviors are in good correspon- 
dence, which is a strong evidence for the validity of our 
assumptions. The experimental relation for N(t) calcu- 
lated in KPZ systems is shown in Fig. || together with 
the curve fitted based on the KAJM result, Eq. (^). 

An associated quantity is the average height of the in- 
terface in time, which can be measured in a very simple 
way in contrast to N{t) mentioned in the preceding para- 
graphs. Let SJt) denote the uncovered fraction of a PNG 
system, Eq. (0). Thus, 



{hKPz{t)) = Y hxpz 
^ Jo 



(x, t) dx 



S{t) dr. (10) 



Fig. H shows the empirical relation fomid for {hxpzit)) 
and the adequately scaled error function for a given in- 
terface evolution. Again, the theory proves to be con- 
sistent with our experimental data to a high degree of 
confidence. 
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FIG. 6. The averaged interface height {h) versus time t, 
the latter measured in lattice updating steps. A system of 
10000 lat tice units is used. The guide line is obtained from 
fitting •^7r/4 j4 erf (_Bf) in a least-squares sense with A = 555 
and B = 1.27 • 10"*. 



In order to verify this assumption, we perform an "is- 
land counting" measurement on the KPZ ensembles with 
the ultimate goal of deriving the number of islands ver- 
sus time statistics, N{t). Islands are interpreted in a 
straightforward way: an island is created in the associ- 
ated PNG system whenever a pinning event occurs, and 
two islands adjoin to form one at the facet summits of 



V. SUMMARY 

In this paper we point out a mapping between in- 
terfaces in two-dimensional systems, described by the 
negative-A quenched Kardar-Parisi-Zhang equation and 
a one-dimensional island growth model. An investigation 
of KPZ interfaces reveals that pinning processes give rise 
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to a morphology with valleys and facets that are formed. 
The facets are found to have straight edges with vanish- 
ing deviations on a macroscopic scale, as well as sharp 
corners. This establishes a connection with the PNG or 
KAJM model. The correspondence is supported by evi- 
dence based on comparing statistical quantities in both 
ensembles. Considering the KPZ model from this point 
of view, it can be stated that pinning events are the mani- 
festation of a local nucleation process that is uncorrelated 
both in time and space. The pinning probability for an 
interface segment which advances covering an area of dA 
only depends on dA in an exponentially decaying fashion. 

The work presented here gives a geometric descrip- 
tion of the A < 0-depinning transition, in analogy to the 
directed percolation depinning-picture of the opposite, 
A > 0-case. There it is possible to extend the correspon- 
dence of the geometrical construction and the pinned in- 
terfaces at Fc to higher dimensions ||2^. In our case this 
still awaits further work, but one should note that the 
KAJM model of island growth becomes much more difh- 
cult to handle, analytically, in d > 1 ]l8| , p^ . 

This work has been supported by the Academy of Fin- 
land's Centre of Excellence Programme. 
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F 


Size 


Pinnings 


Gradient left 


Gradient right 


13.5 


500 
1000 
2500 


10.63 ± 3.02 
10.26 ± 2.07 
9.93 ± 1.24 


1.483 ± 0.087 
1.483 ± 0.066 
1.482 ± 0.039 


1.484 ± 0.086 
1.484 ± 0.067 
1.482 ± 0.040 


14.0 


500 
1000 
2500 


5.80 ± 2.13 
5.65 ± 1.46 
5.53 ± 0.89 


1.652 ±0.085 
1.652 ±0.069 
1.651 ±0.041 


1.654 ±0.085 
1.653 ±0.066 
1.650 ± 0.041 


14.5 


500 
1000 
2500 


3.22 ± 1.42 
2.85 ± 1.02 
2.79 ± 0.63 


1.798 ± 0.093 
1.802 ± 0.068 
1.800 ±0.041 


1.797 ± 0.095 
1.802 ± 0.076 
1.800 ± 0.042 


14.8 


500 
1000 
2500 


2.49 ± 0.94 
1.86 ± 0.80 
1.74 ± 0.50 


1.869 ± 0.085 

1.878 ± 0.073 

1.879 ± 0.039 


1.869 ± 0.093 
1.879 ± 0.069 
1.877 ± 0.039 



TABLE I. Comparison table of interface morphology for close-to-critical systems with different external driving 
forces; Fc ~ 15. The "Pinnings" column contains the average number of valleys found in the systems, normalized 
to 1000 lattice units. The "Gradient left" and "Gradient right" columns show the average slope on the left hand 
side and right hand side of the facets, respectively. The number of interface samples averaged is 10000 in each of 
the cases. The errors are measured by the standard deviation of the respective quantities. 
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